Cavity QED implementation of multi-qubit refined Deutsch-Jozsa algorithm 
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We theoretically study a realization of multi-qubit refined Deutsch-Jozsa (DJ) algorithm using 
resonant interaction of many Rydberg atoms with a single-mode microwave cavity, in which the 
multi-qubit controlled phase gates could be accomplished efficiently. We show how to achieve a 
multi-qubit refined DJ algorithm in high-fidelity even in the case of weak cavity decay. We argue 
that the required operations in our scheme are almost within the present experimental possibilities. 

PACS numbers: 



Quantum algorithms have displayed the unusual power 
in improving computational speed over their classical 
counterparts, due to computational parallelism or inter- 
ference effects. Among the most frequently mentioned 
algorithms the Deutsch-Jozsa (DJ) algorithm is 

the simplest, but demonstrates the power of quantum 
mechanics by distinguishing the constant functions from 
the balanced functions using only one-step logical com- 
putation regardless of the input size. Experimentally, the 
original DJ algorithm Hand/or its modified version (i.e., 
refined DJ algorithm) [J] have been implemented in nu- 
clear magnetic resonance (NMR) system 5, 6], quantum 
dot @, linear optical system jjj], and trapped ions 
Additionally, there have been some theoretical propos- 
als for achieving DJ algorithm using trapped electrons 
[To| . polyatomic molecule s atomic ensembles [l2j], 
Josephson charge qubits |13l| and cavity quantum elec- 
trodynamics (QED) [14l4l( 

In this paper, we will focus on an implementation 
of multi-qubit refined DJ algorithm with the Rydberg 
atoms simultaneously passing through a single-mode mi- 
crowave cavity system. Relevant experiments have been 
carried out for the resonant interaction of two or three 



atoms with cavity mode [17], the entanglement between 
two atoms in a microwave cavity fl8j , and the two-qubit 
phase gate [19(. Based on the available techniques, our 
scheme could have following favorable characters: (i) By 
a smart encoding, we may accomplish the multi-qubit 
gating by one step, which helps to achieve a straight- 
forward and fast implementation of the refined DJ al- 
gorithm. This makes the necessary operations more ef- 
ficient with respect to the previous ideas [H, HU using 
two-qubit conditional gates; (ii) As the DJ algorithm re- 
alized in our scheme is the refined version Q, the auxil- 
iary qubits are not necessary, which might greatly sim- 
plify the experimental requirement particularly in the 
case of scalability; (iii) Although the resonant interac- 
tion is unavoidably affected by the cavity decay, the fast 
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FIG. 1: Quantum circuit for the refined Deutsch-Jozsa algo- 
rithm, where Uf is the /-controlled gate and the operation H 
denotes the Hadamard gate. 



implementation could effectively suppress the detrimen- 
tal effect in the case of the weak cavity decay. 

Let us first briefly review the main idea of the original 
DJ algorithm, which can distinguish constant functions 
/c(x) from balanced functions /s(x) in an A-qubit sys- 
tem using only one query of binary- valued function f{x), 
namely, f(x) : {0,1} N -> {0,1} H. In the original DJ 
algorithm [3(, the function is characterized by the uni- 
tary operation U \x) \y) = \x) \y © f(x)) , where x is an 
A-qubit input but the auxiliary qubit y must be prepared 
in the superposed state (|0) — |l))/\/2, which results in 
the final transformation 



U\x)\y) = U\x)(\0)-\l))/s/2 

->(-l)^)|x)(|0)-|l))/>/2. 



(1) 
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\y) seems superfluous because it keeps unchanged dur- 
ing the operation process, although it actually plays a 
crucial role in above DJ algorithm. This redundancy of 
the auxiliary qubit y is fully removed in the refined DJ 
algorithm [4| , where the action of the /-controlled gate 

U ( f N) can be denoted by U ( f N) \x) = (-1)^ |ar) with A 
being the qubit number. The basic quantum circuit to 
perform this refined DJ algorithm is sketched in Fig. 1. 
For A = 2, operation U { f N) can be specifically expressed 

as (if = diag{{-iy(°°\ (-l)H 01 \ (-1) /(10) , (-1) /(11) } 
in the state space spanned by {|00) , |01) , 1 10} , 1 11)}. It 

" (2) 

has been shown that U) can be reduced to a direct 



2 



product of single-qubit operations if AT < 2. This im- 
plies that the realization of the DJ algorithm could be 
considerably simple. However, for N > 3, the situation 
will become much more intricate. In what follows, we 
will mainly work on a three-qubit implementation of re- 
fined DJ algorithm. Our scheme is directly extendable 
to many-qubit cases. 

In the three-qubit case, along with two /-controlled 
±diag{l, 1, 1, 1, 1, 1, 1, 1} corresponding to 



(a) 



gates UfJ 



the constant functions fc(x), the number of U 
sponding to the balanced functions is C| 



(3) 

fB 



(3) 



Actually, there are only one Uf c 



corre- 
= 70. 
and 35 nontrivial and 



' (3) 

distinct Uyj if we take the symmetry into account and 



neglect the overall phase factors. So our task here is to 
implement 36 unitary transformations. 

As the goal of the DJ algorithm is to differentiate the 
constant functions from the balanced functions, instead 
of finding how Uf N ^ works specifically, we may simply 
consider the case below with one balanced function and 
the corresponding / -controlled operation 



tff^ = diag{l, -1, 1, -1, -1, 1, 1, -1}, 



(2) 



where the state space is spanned by 
{|000) , |001) , |010) , |011) , 1 100) , |101) , |110) , | 111) 
}• 

For clarity of following description, we first consider an 
ideal situation without cavity decay, where three identi- 
cal three-level atoms are input with identical velocities 
and simultaneously interact with the single-mode vac- 
uum cavity field. The atomic internal states are denoted 
by \ij) , \gj) and |e_y) , with decoupled from other 
two states throughout our scheme due to large detuning, 
as shown in Fig. 2 (a). In units of h — 1, the effec- 
tive Hamiltonian in the interaction picture reads Hi = 
7=i ^j(a + 5' J ~ + aS^), where Qj is the coupling con- 

e j) (9j\ 



stant of the jth atom to the cavity mode, Sj~ 



and Sj = \gj) (ej\ are the atomic spin operators for rais- 
ing and lowering, respectively, and a + (a) is the creation 
(annihilation) operator for the cavity mode. 

The effective Hamiltonian Hi is similar to the quan- 
tum computing model of trapped ions in linear trap 
poj . where the multi-qubit controlled phase flip (CPF) 
gate J p = I — 2 \p) (p\ with / the identity matrix and 
\p) = |000) , |001) , • • •, and |111) , could be achieved with 
a high success probability and a high fidelity In the 
present scheme, we will try to move such an idea to cavity 
QED system for designing CPF gate, which is essential 
to the / -controlled operation Uj . Different from the 
time- varying coupling strength f2,-(i) in (20l |. however, the 
coupling strength in our scheme is time- independent, i.e., 
0, = fio cos(27T2/Ao) exp(— r 2 /w 2 ) ~ f2 cos(27rz/Ao) 
[2 If , where r is the distance of the atom away from 
the cavity center, and Ao and w are the wavelength 
and the waist of the cavity mode,respectively. We use 
the same qubit definitions as in [2fJ, that is, the logic 
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FIG. 2: (a) Schematic setup for implementing three-qubit 
refined DJ algorithm in the case of the /-controlled gate is 
Uf^ = diag{l, — 1, 1, — 1, — 1, 1, 1, —1}, where the inset shows 
the atomic level structure, and Di, D2, D3 are state-selective 
field-ionization detectors, (b) The detailed configuration of 
the unitary operations J ei i 2 i 3 , Je ig2 g 3 > Jg\i 2 i 3 an d ^9192^3- 
The cavity is a microwave cavity sustaining a single mode 
with a standing-wave pattern along the 2-axis. The atoms 1, 
2 and 3 prepared in high-lying circular Rydberg states are sent 
through the cavity with proper speed, resonantly interacting 
with the cavity mode. Ri(R i ) with i = 1, 2, 3 denote Ramsey 
zones for performing single-qubit rotations as explained in the 
text. 



state |0) (|1)) of the qubit 1 is denoted by \g{) (|ei)) 
of the atom 1; |<7a) and \i 2 ) of the atom 2 encode the 
logic state |0) (|1)) of the qubit 2; the logic state |0) 
(|1)) of the qubit 3 is represented by \g 3 ) (|ia)) of the 
atom 3. If we assume the different atoms with coupling 
constants fii : Cl 2 : O3 = 1 : 10 : 10 and the gating time 
To = 7r/f2i, we can obtain an approximate three-qubit 
CPF gate Jin = J ei i 2 i 3 = diag{l, 1, 1, 1, at\, a 2 , a 3 , -1} 
in the computational subspace spanned by 

\gi) I32) I33) , \gi) Ifla) 1*3) , \gi) \h) \ga) , \gi) V2) \h) , 
|ei) |fl2> \9a) , |ei) \g 2 ) fa) , |ei) \i 2 ) \ga) , |ei) |i 2 ) \h) , 
where a; = [f^cos(9,) + 9? - fi?]/9f 1 with 
i = 1,2,3, and 9i = y^f 

9 3 = y/sti + n 2 3 . 



nl 9 2 = 



Based on the gate J ei i 2 i 3 , we construct the /-controlled 

~ (3) 

gate Uf by a straightforward way, i.e., 



where J e 



r}(3) — j j 7 7 

u f B i ~ Je i l 2»3 u e 1 g 2 g 3 J g\i 2 i 3 ° g\g 2 i 3 i 
'253 Jl00 &x,3&x,2Jeii 2 i 3 &x,2&x,3j ^ g\i 2 i 3 



(3) 
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Jon — S Xt iJ ei i 2 i 3 S Xt i and Jg x g 2 i3 — Jom — 
Vx,2S Xt iJ eii2 i 3 S Xtl a Xt 2,vfitYi cr xJ = \ij){9j\ + \gj){ij\ with 
j 7^ 1, and S x .\ = |ei)(gi| + |<?i)(ei|. The other four indis- 
pensable CPF gates could also be created as J ei i 2 g 3 — 

JlXO ^x^J e.\i 2 %z@x,3i Jeyg-ziz — <-^101 — ®x$J&\iiis®x % %i 
Jgii293 = ^010 = ^x,3S Xl lJ e ,i 2 i 3 S Xt iax,3i Jgig 2 93 = ^000 = 
&X,3&X,2$ Xi lJ 61*2^3 ^X z i&x,2&x,3- Using these eight CPF 

gates, we could construct other 35 /-controlled gates 



U 



(3) 



(m = 1, 2, • • •, 35), with each U^J involving four 
different CPF gates. 

Along with the above-mentioned CPF gates, three- 
qubit Hadamard gates H® 3 = n*=i H should be per- 
formed to encode the input and decode the output, re- 
spectively, as depicted in Fig. 2(a), where Hi — [(\0) i + 

1 1 > i ) (0 1 i + ( 1 0) t - 1 1 ) { ) ( 1 1 i } j V2 is the Hadamard gate act- 
ing on the ith atom. These gates could be achieved using 
external microwave pulses. As a result, a full three-qubit 
refined D J algorithm is available. Taking the /-controlled 
gate Uj gi as an example, we have designed a three-qubit 
refined DJ algorithm setup in Fig. 2, where three Ry- 
dberg atoms prepared in the state |\&o) = \g1g2g3) are 
initially encoded by the three-qubit Hadamard gate H® 3 
[23| . and then sent through the cavity with the identi- 
cal velocities. The implementation could be simply de- 
scribed as, 

|*o) H® 3 U { ^ m H m 



|* /(n ) = l{A m |001) + B m |010) + C m |011) 
+ £ m |101) + F m |110) + G m |lll)}; 



47) = ± 1 000) , 



Dm 1 100) 

(4) 

(5) 



For an initial state \*S?q) = |000), after the opera- 
tions in Fig. 2, the output state of the three atoms 



in the decay case is given by 

A m |001) +B m |010) +C m |011> - 
|110) + G m |lll), where R m p 
slightly deviated 
demonstrate the 
tionF m = W fm m' 



= Rn 



r fm/ - ,. m |000) + 

^ m |100)+i? m |101) + 
F m 1 110) + G m I 111), where R m w and A' m ,B-' m , ■■■,G m 
are slightly deviated from A m , B mi ■ ■ -,G m . In Fig. 3(a), 
we demonstrate the fidelity according to following rela- 

/ /™l*/m)(*/ml*/'») 25]' in which we have 
also considered the influence from the deviation A due 
to slightly different atomic velocities. 

In addition to the imperfection considered above, there 
are other noise effects need to investigate, such as the 
resonant dipole interaction between two neighboring Ry- 
dberg atoms. Since the typical dipole moment of Ryd- 
berg atoms is about several hundred q e ao, where ao is 
the Bohr radius and q e is electron charge, we can make 
a calculation of the dipole coupling strength S between 
two neighboring Rydberg atoms being in the same range 
as the cavity decay rate ~ KHz. Considering the system 
Hamiltonian (Eq. (6)), we should also take an additional 

term J2j=i ^\ e j e j+i)( e j e j+i\ m t° account to assess the 
influence of the nearest-neighbor dipole-dipole interac- 
tion. 

In the present of cavity decay, we have plotted in Fig. 

" (3) 

4 the fidelity of Uf^ in a three-qubit refined DJ algo- 
rithm versus the parameter 6/Q1 and different cavity de- 
cay rates. As shown in Fig. 4, as long as these dipole 
interactions are weak, our scheme can still achieve a high 
fidelity. In a realistic experiment, the situation would be 
more complicated than our consideration above. So to 
carry out our scheme efficiently and with high fidelity, we 
have to suppress these above-mentioned imperfect factors 
as much as we can. 



where £ is a normalized coefficient with A m , B m , • • •, and 
G m being or ±1. It implies that, if f(x) is constant, 
the state of the atoms becomes |000) ; But if f(x) is bal- 
anced, the state of the atoms becomes a superposition 
state \^>f) , excluding the component |000) . So we can 
efficiently determine whether the function is constant or 
balanced by a collective measurement on the output state 
of the three atoms. 

Under the assumption of weak cavity decay that no 
photon actually leaks out of the microwave cavity during 
our implementation, we may reconsider our scheme using 
the quantum trajectory method [2l|, 



Hr 



E 



^(a+Sr+aSf)-^ 



(6) 



where k is the cavity decay rate. When we choose the 
atom-cavity interaction time Tjj = tt/Ai k with Ai K = 
yOf — k 2 /16 and meet the condition 57i : Q2 : ^3 = 1 : 
10 : 10, the approximate three-qubit CPF gate in the 
decay case becomes J ei j 2 j 3 = diag{l, 1, 1, 1, c±, ci, C3, C4} 
[24| , and the gating time becomes Tp = it/ ^JVt\ — k 2 /16. 




FIG. 3: (a) Fidelity of (jf^ in a three-qubit refined DJ algo- 
rithm versus the parameter /t/fii and the deviation of gating 
time A. (b) Fidelity of llf^ in a four-qubit refined DJ algo- 
rithm versus the coupling ratio 77 = fii/fi and the parameter 



4 



Extending to the many-qubit case, we could construct 

the iV-qubit CPF gate Vj m i) = diag{l 7 1, 1, • • •, -1} 

by meeting the condition for coupling constants 
Qi <C r?2 = ^3 = • • • = = ^ and by keeping the 
gating time To = n/ y/flf — n 2 /16 unchanged. This 
implies that the CPF gate could be carried out by a 
constant time irrelevant to the qubit number, which is 
favorable for a scalable DJ algorithm in cavity QED 
system. Since the single-qubit operation takes negligible 
time in comparison with that for the many-qubit CPF 
gate, we may roughly omit the single-qubit operation 
time, which yields running time of the A-qubit refined 
DJ algorithm to be T {N) = 2 N ~ 1 x T D . Fig. 3(6) demon- 
strates an example of four-qubit /-controlled gate uj^ = 
diag{l, -1,1, -1,-1, 1,-1, 1,1, -1,-1, 1,1, -1,-1,1}, 
constructed by 




0.992 L ■ ■ ■ ■ J 

0.00 0.02 0.04 0.06 0.08 0.10 

sia. 



" (3) 

FIG. 4: Fidelity of U)^ in a three-qubit refined DJ algorithm 
versus the parameter 6/fli and the cavity decay rate k, where 
the solid, dotted, and dashed curves represent the cases of 
k = ili/10, K = fli/20, and k = fh/50, respectively. 

We briefly address the experimental feasibility of our 
scheme by considering three high-lying Rydberg atomic 
levels with principal quantum numbers 49, 50 and 51 to 
be levels \g), \i) and |e) , respectively. Based on the exper- 
imental numbers reported in [2o| . the coupling strength 
at the cavity centre could be fig = 27r x 51 KHz, and 
the Rydberg atomic lifetime is 30 ms. Specifically, as- 
suming Q = 2n x 51 KHz, fli = 2n x 5.1 KHz, and 
k = 10~ 3 £li, direct calculation shows t hat the time for a 
single A-qubit CPF gate is T D = Tr/s/T% - k 2 /16 w 98 
/is, which is much shorter than either the cavity decay 
time, i.e., 2-k/k s» 0.2 s, or the Rydberg atomic lifetime. 
Actually, the lifetime of the photon in the superconduct- 
ing cavity has reached 0.5 s recently [26]. To make the 
quantum trajectory menthod workable, we require our 
implementation time to be much shorter than the cav- 
ity decay time, which yields the condition k ^ fli/2 N ~ 2 . 
Considering the values listed above, we find that N could 
be 9. 

In current microwave cavity experiments (2lj . the 
intra-atom interaction occurs in the central region of the 



cavity with the Rabi frequency Slj — Hq cos(2-7rz/Ao) 
To meet the condition of Qi = fl/10, the A 
atoms should be sent through the cavity with the first 
atom going along the y-axis deviating from the nodes 
by (Ao/27r) arccos(l/10), but other atoms through the 
antinodes. Specifically, for the A input atoms (suppose 
iV being an even number for simplicity), we have the 
tracks of the atomic movement as z\ = — (A/2)Ao + 
(A /27r)arccos(l/10), z 2 = -(N/2 - 1)A , z N = 
(N/2)\q. This is of course highly challenging with cur- 
rent experimental technology because we have not yet 
found any experimental report for A 3 atoms si- 
multaneously going through a microwave cavity. How- 
ever, the two-atom entanglement in a microwave cavity 
has been achieved using van der Waals collision between 
the atoms [18j under a non-resonant condition. More- 
over, our scheme could also be straightforwardly applied 
to other quantum information processing candidate sys- 
tems, e.g., the ion-trap-cavity combinatory setup 12711 or 
cavity-embedded optical lattices confining atoms [28[, in 
which the atoms are localized very well. 

Alternatively, the superconducting circuit QED [29l — 
[3lT ] would become a more suitable candidate for imple- 
menting our scheme. In the circuit QED [2!|[3(|, a num- 
ber of superconducting qubits play the role of artificial 
atoms and quantum bus is provided by a ID supercon- 
ducting transmission line resonator. So the dynamics of 
a two-level qubit coupled to a single mode of an electro- 
magnetic field is also described by the Jaynes-Cummings 
Hamiltonian, which is an essential requirement for the 
present scheme. Implementing multi-qubit refined DJ 
algorithm in the circuit QED has several merits that are 
worth mentioning here. First, the experiments performed 
in the circuit QED set-up have demonstrated that long 
coherence time and very strong coupling could be real- 
ized [3(| , which could greatly reduce the gating time and 
noise effects in our scheme. Second, the widely separated 
superconducting qubits are well localized in the circuit 
QED, so full in-situ control over the qubit parameters 
and system qubit Hamiltonian is easily achieved. Third, 
the strong coupling between the field in the resonator 
and the qubits can be used to perform a single-shot high 
efficiency quantum nondemolition (QND) readout of the 
state of the qubits without the need for additional signal 
ports. 

In conclusion, we have proposed a potentially practical 
scheme for realizing a multi-qubit refined DJ algorithm 
by resonant interaction of Rydberg atoms in a microwave 
cavity. We have estimated the influence from the cavity 
decay on our scheme and shown that our scheme could be 
achieved efficiently to distinguish the balanced functions 
from the constant functions with high fidelity. We argue 
that our present scheme would be helpful for demonstra- 
tion of refined DJ algorithm at large scale using cavity 
QED devices. 

This work is supported by NNSF of China under No. 
10774163 and No. 10774042. 
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